Controllable coherent population transfers in superconducting qubits for quantum computing 



oo 

o 

O 

(N 



Oh, 



> 

o 

00 
O 



X 



L.F. Wei/' 2 J.R. Johansson,^ L.X. Cen/ S. Ashhab,^-^ and Franco Nori'-^-^ 

'crest, Japan Science and Technology Agency (JST), Kawaguchi, Saitama 332-0012, Japan 
^Laboratory of Quantum Opt-electronics, Southwest Jiaotong University, Chengdu 610031, China 
Frontier Research System, The Institute of Physical and Chemical Research (RIKEN), Wako-shi, Saitama, 351-0198, Japan 
'* Department of Physics, Sichuan University, Chengdu, 610064, China 

^Center for Theoretical Physics, Physics Department, CSCS, 
The University of Michigan, Ann Arbor, Michigan 48109-1040, USA 
(Dated: February 3, 2008) 

We propose an approach to coherently transfer populations between selected quantum states in one- and 
two-qubit systems by using controllable Stark-chirped rapid adiabatic passages (SCRAPs). These evolution- 
time insensitive transfers, assisted by easily implementable single-qubit phase-shift operations, could serve as 
elementary logic gates for quantum computing. Specifically, this proposal could be conveniently demonstrated 
with existing Josephson phase qubits. Our proposal can find an immediate application in the readout of these 
qubits. Indeed, the broken parity symmetries of the bound states in these artificial "atoms" provide an efficient 
approach to design the required adiabatic pulses. 



PACS number(s): 42.50.Hz, 03.67.Lx, 85.25.Cp. 

Introduction. — The field of quantum computing is attract- 
ing considerable experimental and theoretical attention. Usu- 
ally, elementary logic gates in quantum computing networks 
are implemented using precisely designed resonant pulses. 
The various fluctuations and operational imperfections that 
exist in practice (e.g., the intensities of the applied pulses 
and decoherence of the systems), however, limit these de- 
signs. For example, the usual 7r-pulse driving for performing 
a single-qubit NOT gate requires both a resonance condition 
and also a precise value of the pulse area. Also, the difficulty 
of switching on/off interbit couplings |1] strongly limits the 
precise design of the required pulses for two-qubit gates. 

Here we propose an approach to coherently transfer the 
populations of qubit states by using Stark-chirped rapid adi- 
abatic passages (SCRAPs) As in the case of geometric 
phases |3], these population transfers are insensitive to the 
dynamical evolution times of the qubits, as long as they are 
adiabatic. Thus, here it is not necessary to design before- 
hand the exact durations of the applied pulses for these trans- 
fers. This is a convenient feature that could reduce the sen- 
sitivity of the gate fidelities to certain types of fluctuations. 
Another convenient feature of our proposal is that the phase 
factors related to the transfer durations (which are important 
for the operation of quantum gates) need only be known after 
the population transfer is completed, at which time they can 
be cancelled using easily implementable single-qubit phase- 
shift operations. Therefore, depending on the nature of fluc- 
tuations in the system, rapid adiabatic passages (RAPs) of 
populations could offer an attractive approach to implement- 
ing high-fidelity single-qubit NOT operations and two-qubit 
SWAP gates for quantum computing. Also, the SCRAP-based 
quantum computation proposed here is insensitive to the ge- 
ometric properties of the adiabatic passage paths. Thus, our 
approach for quantum computing is distinctly different from 
both adiabatic quantum computation (where the system is al- 
ways kept in its ground state [4]) and holonomic quantum 
computating (where implementations of quantum gates are 



strongly related to the topological features of either adiabatic 
or non-adiabatic evolution paths Jst]). 

Although other adiabatic passage (AP) techniques, such as 
stimulated Raman APs (STIRAPs) [6], have already been pro- 
posed to implement quantum gates ||7[], the present SCRAP- 
based approach possesses certain advantages, such as: (i) it 
advantageously utilizes dynamical Stark shifts induced by the 
applied strong pulses (required to enforce adiabatic evolu- 
tions) to produce the required detuning-chirps of the qubits, 
while in STIRAP these shifts are unwanted and thus have to 
be overcome for performing robust resonant drivings; and (ii) 
it couples qubit levels directly via either one-photon or mul- 
tiphoton transitions, while in the STIRAP approach auxiliary 
levels are required. 

The key of SCRAP is how to produce time-dependent 
detunings by chirping the qubit levels. For most natu- 
ral atomic/molecular systems, where each bound state pos- 
sesses a definite parity, the required detuning chirps could be 
achieved by making use of the Stark effect (via either real, but 
relatively-weak, two-photon excitations of the qubit levels_|8|] 
or certain virtual excitations to auxiliary bosonic modes [9]). 
Here we show that the breaking of parity symmetries in the 
bound states in current-biased Josephson junctions (CBJJs) 
provides an advantage, because the desirable detuning chirps 
can be produced by single-photon pulses. This is because all 
the electric-dipole matrix elements could be nonzero in such 
artificial "atoms" [10]. As a consequence, the SCRAP-based 
quantum gates proposed here could be conveniently demon- 
strated with driven Josephson phase qubits 111 111 generated by 
CBJJs. In order to stress the analogy with atomic systems, 
we will refer to the energy shifts of the CBJJ energy levels 
generated by external pulses as Stark shifts. 

Models. — Usually, single-qubit gates are implemented by 
using coherent Rabi oscillations. The Hamiltonian of such a 
driven qubit reads i/o(i) = t^oCz/S + R{t)ax, with ujq being 
the eigenfrequency of the qubit and R{t) the controllable cou- 
pling between the qubit states; and ax are Pauli operators. 
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If the qubit is driven resonantly, e.g., R{t) = ^{t) cos{u!Qt), 
then the qubit undergoes a rotation Rx{t) = cos[A{t)/2] — 
ia^ sm[A{t)/ 2], v/ith A(t) = J*n(t')dt'. For realizing a 
single-qubit NOT-gate, the pulse area is required to be pre- 
cisely designed as A{t) = tt, since the population of the 
target logic state P{t) = [1 — cosA{t)]/2 is very sensitive 
to the pulse area A{t) [in this example, we are assuming an 
initially empty target state]. Relaxing such a rigorous con- 
dition, we additionally chirp the qubit' s eigenfrequency 
by introducing a time-dependent Stark shift A(t). There- 
fore, the qubit evolves under the time-dependent Hamiltonian 
iJg(i) = uioaz/2 + R{t)<Tx + A{t)(Tz/2, which becomes 



Hi{t) 



1/0 n{t) \ 

2 [ n{t) 2A{t) ) 



(1) 



in the interaction picture. Under the condition 

<^[A\t) + n\t)]^/^, (2) 



dt 



dt 



the driven qubit adiabatically evolves along two paths — the in- 
stantaneous eigenstates |A_(t)} = cos[6'(t)] |0) — sin[0(i)]|l) 
and \\+{t)) = sin[e'(i)]|0) + cofi[0{t)]\l), respectively. In 
principle, these adiabatic evolutions could produce arbitrary 
single-qubit gates. For example, a single detuning pulse 
A{t) (without a Rabi pulse) is sufficient to produce a phase- 
shift gate: Uz{a) = cxp(ja|l)(l|), a = - J^^ A{t)dt. 
Furthermore, combining the Rabi and detuning pulses for 
rotating the mixing angle d{t) — arctan[il(t)/A(t)]/2, 
from 9{—oo) = to d{+(x) — tt/2, another single- 
qubit gate Ux — exp(i/3+)CT+ — exp{i(3_)a^ (with (3± = 

~ A^±(*)rf^> M±(^) = A(t) ± ^A^t) + n^t)) can be 
adiabatically implemented as: 



Ux 



|A_(-(X3)) = |0) 



|A-(t)) 



|A+(-oo)) = |1) |A+(+oo)) = e'/5+|0). 



(3) 



This is a single-qubit rotation that completely inverts the pop- 
ulations of the qubit's logic states and thus is equivalent to the 
single-qubit NOT gate. Note that here the population trans- 
fer is insensitive to the pulse duration and other details of the 
pulse shape — there is no need to precisely design these be- 
forehand. Different durations for finishing these transfers only 
induce different additional phases f3±, which can then be can- 
celled by properly applying the phase shift operations Uz{ct). 

Similarly, the applied pulses are usually required to be 
exactly designed for implementing two-qubit gates. For 
example il2[l . for a typical two-qubit system described 
by the XY-type Hamiltonian H12 = X]i=i 2 "^j^^z^V^ + 

K{t)J2i=£j=i2'^'+'^'^^/'^^ with switchable real interbit- 
coupling coefficient K{t), the implementation of a two-qubit 
SWAP gate requires that the interbit interaction time t should 
be precisely set as /g K{t')dt' = vr (when uji = cj2)- This 
difficulty could be overcome by introducing a time-dependent 
dc-driving to chirp the levels of one qubit. In fact, we can add 



ID 






cn 


|1> 










I1> 


CD 
tz 
LLI 


|0> 




Z' ii> 
'l^-(t)> 




t 








t 





FIG. 1: (Color online) Simulated SCRAPs for inverting the qubit's 
logic states by certain pulse combinations: (left) a linear detuning 
pulse A(t) = VaX, combined with a constant Rabi pulse f2(t) = Sl^; 
and (right) a linear detuning pulse A(t) = vi,t, assisted by a 
Gaussian-shape Rabi pulse Q,{t) — fib cxp{—t'^ /T^). Here, the 
solid (black) lines are the expected adiabatic passage paths, and the 
dashed (red) lines represent the unwanted Landau-Zener tunnelling 
paths. 



(2) 

a Stark-shift term A2 (i)o'i /2 applied to the second qubit and 
evolve the system via (in the interaction picture) 



( -A2(t) 
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K(t) 





K{t) 
A2(i) 

A2(i) ) 
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Obviously, three invariant subspaces; 5Ro — {1 00)}, 5Ri — 
{|11)}, and K2 = {|01), |10)} exist in the above driven dy- 
namics. This implies that the populations of states |00) and 
1 11) are always unchanged, while the evolution within the 
subspace 5R2 is determined by the reduced time-dependent 
Hamiltonian (1) with ri(<) and A(t) being replaced by K(t) 
and A^ii), respectively. Therefore, the APs determined by 
the Hamiltonian H.2{t) produce an efficient two-qubit SWAP 
gate; the populations of |00) and |11) remain unchanged, 
while the populations of state |10) and |01) are exchanged. 
The passages are just required to be adiabatic and again are 
insensitive to the exact details of the applied pulses. 

Figure 1 shows schematic diagrams of two simulated 
SCRAPs. There, solid (black) lines are the desirable AP paths, 
and the dashed (red) lines are the unwanted Landau-Zener tun- 
nellings 1 13] (whose probabilities should be negligible for the 
present adiabatic manipulations). These designs could be sim- 
ilarly used to adiabatically invert the populations of states 1 10) 
and |01) for implementing the two-qubit SWAP gate. 

Demonstrations with driven Josephson phase qubits. — In 
principle, the above generic proposal could be experimentally 
demonstrated with various physical systems e.g., the gas- 
phase atoms and molecules, where SCRAPs are experimen- 
tally feasible. Here, we propose a convenient demonstration 
with solid-sXsit Josephson junctions. 

A CBJJ (see, e.g., [11]) biased by a time-independent dc- 
current /{, is described by i/o — /2m + U{Ib,6). For- 
mally, such a CBJJ could be regarded as an artificial "atom", 
with an effective mass m = Cj$o/(27r), moving in a po- 
tential [/(/(,, (5) ~ ~Ej{cos6 — IbS/Io). Here, Iq and 
Ej = ^qIo/2tt are, respectively, the critical current and the 
Josephson energy of the junction of capacitance Cj. Under 
proper dc-bias, e.g., /(, < Iq, the CBJJ has only a few bound 
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FIG. 2: (Color online) SCRAP-based population transfers in a sin- 
gle Josephson phase qubit. (Left) manipulated scheme: CBJJ levels 
with dashed chirped qubit energy splitting A{t) is coupled (solid ar- 
row) by a Rabi pulse ^l{t). Dotted red arrow shows the unwanted 
leakage transition between the chirping levels |1) and j2). (Right) 
time-evolutions Pj(t) of the occupation probabilities of the lowest 
three levels |j> (j = 0, 1, 2) in a CBJJ during the designed SCRAPs 
for inverting the populations of the qubit logic states. This shows that 
during the desirable SCIiAPs the qubit leakage is negligible. 



states: the lowest two levels, |0) and |1), encode the qubit of 
eigenfrequency wio = {Ei — Eo)/h. During the manipula- 
tions of the qubit, the third bound state |2) of energy E2 might 
be involved, as the difference between E2 — Ei and Ei — Eq is 
relatively small. Due to the broken mirror symmetry of the po- 
tential U{6) for S ^ —S, bound states of this artificial "atom" 
lose their well-defined parities. As a consequence, all the 
electric-dipole matrix elements 5ij = = 0,1,2, 

could be nonzero [10]. This is essentially different from 
the situations in most natural atoms/molecules, where all the 
bound states have well-defined parities and the electric-dipole 
selection mle forbids transitions between states with the same 
parity. By making use of this property. Fig. 2 shows how 
to perform the expected SCRAP with a single CBJJ by only 
applying an amplitude-controlled dc-pulse /dc(i) (to slowly 
chirp the qubit's transition frequency) and a microwave pulse 
Iac{t) = Aoi{t) cos{ujoit) (to couple the qubit states). Under 
these two pulses, the Hamiltonian of the driven CBJJ reads 
Hi{t) ^ Ho- ($o/27r)[/dc(t) + Iac{t)]5. Neglecting leak- 
age, we then get the desirable Hamiltonian (1) with A(t) = 
A(i) = -($o/27r)/rf,(t)(5n - ^oo) and ^l{t) = ^l{t) - 
— (<l>o/27r)74oi(i)(5oi. Obviously, for a natural atom/molecule 
with 5ii = 0, the present scheme for producing a Stark shift 
cannot be applied. 

Specifically, for typical experimental parameters [fTll] 
(Cj = 4.3 pF, /o 13.3 AiA and 4 = 0.9725/o), our numer- 
ical calculations show that the energy-splittings of the lowest 
three bound states in this CBJJ to be = 5.981 GHz and 
LO21 = 5.637 GHz. The electric-dipole matrix elements be- 
tween these states are Jqo = 1.406, Sn = 1.425, 822 = 
1.450, (5oi = (5io 0.053, 812 = 821 = 0.077, and ^02 = 
820 = —0.004. If the applied dc-pulse is a linear function 
of time (i.e., /dc(i) — vit with vi constant) and the coupling 
Rabi amplitude n{t) = fii is fixed, then the above SCRAP re- 
duces to the standard Landau-Zener problem [13]. For a typi- 
cal driving with vi = 0.15 nA/ns and Aqi = 1.25 nA, Fig. 2 
simulates the time evolutions of the populations in this three- 
level system during the designed SCRAPs. The unwanted (but 
practically unavoidable) near-resonant transition between the 
chirping levels |1) and |2} (due to the small difference be- 



tween llJ21 and uiq) has been considered. Figure 2 shows that 
during the above passages the leakage to the third state |2} is 
sufficiently small. Thus, the above proposal of performing the 
desirable SCRAPs to implement single-qubit gates should be 
experimentally robust. 

The adiabatic manipulations proposed above could also 
be utilized to read out the qubits. In the usual readout ap- 
proach fn\, the potential barrier is lowered fast to enhance 
the tunneling and subsequent detection of the logic state |1). 
Recently UJl, a TT-puIse resonant with the 1 1) ^ |2) transition 
was added to the readout sequence for improved fidelity; The 
tunnelling rate of the state 1 2) is significantly higher than those 
of the qubit levels, and thus could be easily detected. The 
readout scheme used in (l^ can be improved further by uti- 
lizing the above SCRAP by combining the applied microwave 
pulse and the bias-current ramp. The population of state |1) 
is then transferred to state |2) with very high fidelity. In con- 
trast to the above APs for quantum logic operations, here the 
population transfer for readout is not bidirectional, as the pop- 
ulation of the target state |2) is initially empty. The fidelity of 
such a readout could be very high, as long as the relevant AP 
is sufficiently fast compared to the qubit decoherence time. 

Similarly, SCRAPs could also be used to implement two- 
qubit gates in Josephson phase qubits. With no loss of 



generality, we consider a superconducting circuit 111 111 pro- 
duced by capacitively coupling two identical CBJJs. The 
SWAP gate is typically performed by requiring that the two 
CBJJs be biased identically (yielding the same level struc- 
tures) and the static interbit coupling between them reach 
the maximal value K^. If one waits precisely for an in- 
teraction time T — tt/2K(), then a two-qubit SWAP gate 
is produced ilTsll . In order to relax such exact constraints 
for the coupling procedure, we propose adding a control- 
lable dc current, I^J (t) = V2t, applied to the second CBJJ. 
Thus one can drive the circuit under Hamiltonian Hi2{t) — 

Efc=i,2 Hok + {'2n/^ofpip2/C,n-{^a/27T)I^l\t)82. Here, 
the last term is the driving of the circuit, and the first term 
Hok = {27T/^o)^pl/i2Cj)-Ejcos8k~i^o/2Tr)h8kisthe 
Hamiltonian of the fcth CBJJ with a renormalized junction- 
capacitance Cj = Cj(l + C), with C — Cm/{Cj + Cm). 
The coupling between these two CBJJs is described by the 
second term with = C,/[C,j{l + Q] being the effec- 

tive coupling capacitance. Suppose that the applied driving 
is not too strong, such that the dynamics of each CBJJ is still 
safely limited within the subspace 0fc = {|Ofe}, jlfc), |2fc)}: 
Y^i=o\^k){h\ = 1- The circuit consequently evolves 
within the total Hilbert space = 0i 02. Using the 
interaction picture defined by the unitary operator Uq — 
nA;=i 2 exp(— Ym=o l^fc) (^fe|)' we can easily check that, for 
the dynamics of the present circuit, three invariant subspaces 
(relating to the computational basis) exist: (i) 3i = {|00)} 
corresponding to the sub-Hamiltonian Hi — £'oo(i)|00)(00| 

with£;ooW = -[$o/(27r)]4^'(t),5oo + (27r/$o)'pgo/^m' 
Pw = (h\Pk\l'k) and 8u' - {lk\8k\l'k); (ii) ^2 = {|01), |10)} 
corresponding to the sub-Hamiltonian H2 (t) taking the form 
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FIG. 3: (Color online) SCRAPs within the invariant subspace 
Sfa = {102), |11), 1 20)} for the dynamics of two identical three-level 
capacitively-coupled CBJJs driven by an amplitude-controllable dc- 
pulse. (Left) Adiabatic energies and the desirable AP path (the mid- 
dle solid-line with arrows): A — » Ci — » C2 — > C3 B. (Right) 
Time evolutions of populations Pait), ct — 20, 11,02, within the 
invariant subspace ^3 during the designed SCRAPs for inverting the 
populations of |10) and |01). It is shown that the initial population 
of the 1 11) state (corresponding to the A-regime in the left figure) 
is adiabatically partly transferred to the two states |20) and |02) in 
the Ci, C2, and C3 regimes, respectively. Note that the population 
of the state |11) vanishes att — and completely returns after the 
passages. 

of Eq. (1) with n{t) = 2(27r/$o)^P?o/C'm and A{t) 
A{t) = [0o/(27r)]/£^(i)(<5ii - <5oo); and (iii) ^3 = {|02) = 
a), |11) — \b), |20) — \c)} corresponding to 

/ Ea{t) nabC-^'^ \ 

with Ea{t) = -[$o/(2^)]4c^(i)'522 + (2^/$o) 

Eb{t) = -[$o/(27r)]4^^(t)<5ii + (2V$o)Vi/C™,£;cW = 

-[$o/(27r)]4^^(t)5oo + (27r/$o)Woo/C™; ^ab = f^ba = 

flbc = i^cb = (27r/$o)^POlPl2/C'™, ^ac = ^ca = 

(27r/$o)^Po2/C'm, and ■& = wio — ^21- Under the APs 
for exchanging the populations of the states |10) and |01), 
we can easily see that the population of |00) remains un- 
changed. Also, after the desired APs, the population of the 
state 1 11) should also be unchanged. Indeed, this is verified 
numerically in Fig. 3 for the typical parameters ( — 0.05 and 
V2 = 3.0 nA/ns. Therefore, the desirable two-qubit SWAP 
gate could also be effectively produced by utilizing the pro- 
posed SCRAPS. 

Discussions and Conclusions. — By using SCRAPs, we 
have shown that populations could be controllably transferred 
between selected quantum states, insensitive to the details 
of the applied adiabatic pulses. Assisted by readily imple- 
mentable single-qubit phase shift operations, these adiabatic 
population transfers could be used to generate universal logic 
gates for quantum computing. Experimentally existing su- 
perconducting circuits were treated as a specific example to 
demonstrate the proposed approach. 

Like other RAPs, the adiabatic nature of the present 
SCRAPs requires that the passages should be sufficiently slow 
(compared to the usual Rabi oscillations) and sufficiently fast 
(compared to the decoherence times of the qubits). Satisfy- 



ing both conditions simultaneously does not pose any seri- 
ous difficulty with typical experimental parameters. Indeed, 
as we have shown above, experimentally feasible APs could 
be applied within tens of nanoseconds. This time interval 
is significantly longer than the typical period of an experi- 
mental Rabi oscillation, which usually does not exceed a few 
nanoseconds, and could be obviously shorter than the typi- 
cal decoherence times of existing qubits, which might reach 
hundreds of nanoseconds, e.g., for the Josephson phase qubits 
reported in 1 1 1]. Solid-state qubits offer evident advantages 
due to their scalability and controllability. Therefore, RAPs 
in solid-state qubits could provide an attractive approach for 
data storage and quantum information processing. We hope 
that such techniques will be experimentally implemented in 
the near future. 
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